Atomic-phase interference devices based on ring-shaped Bose-Einstein condensates: 

Two ring case 
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We theoretically investigate the ground-state properties and quantum dynamics of a pair of ad- 
jacent ring-shaped Bose-Einstein condensates that are coupled via tunneling. This device, which is 
the analogue of a symmetric superconducting quantum interference device, is the simplest version 
of what we term an Atomic-Phase Interference Device (APHID). The two-ring APHID is shown to 
be sensitive to rotation. 
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INTRODUCTION 



The last few years have witnessed magnificent ad- 
vances in the preparation, manipulation, and exploration 
of atomic Bose-Einstein condensates (BECs). These 
quantum-degenerate systems offer an excellent experi- 
mental platform from which to study a multitude of non- 
linear matter-wave phenomena including four- wave mix- 
ing [Q, dark |Q and bright [|| solitons, superfluid vortices 
Q, and the generation and study of quantized vortices 
on toroidal atomic traps or rings. In particular, ring- 
shaped BECs allow for the study of phenomena related 
to persistent currents and rotational motion, with po- 
tential applications to rotation sensing. In this paper, 
our goal is to take the first theoretical steps in study- 
ing Josephson coupling between adjacent ring BECs (as 
opposed to concentric ring BECs that have been consid- 
ered previously [^). In particular, we investigate how 
quantum tunneling between two condensates trapped in 
adjacent toroidal traps, formed for example using optical- 
dipole traps with Laguerre-Gaussian light beams, mod- 
ifies both the ground state properties and quantum dy- 
namics of the system. The two-ring BEC system is the 
simplest example of what we refer to as an Atomic-Phase 
Interference Device (APHID), essentially a neutral- atom 
analog of a SQUID (Superconducting Quantum Interfer- 
ence Device) . The properties of the APHID will be shown 
to be strongly influenced by the individual phases of the 
matter- waves in the rings. 



The remainder of this paper is organized as follows: In 
the next section we elucidate the details of the model 
we use. Following this, we explore the properties of 
the ground state and first excited state of the system. 
We then look at the Josephson coupling and the time- 
dependent solutions, highlighting important considera- 
tions due to the effects of rotation, followed by concluding 
remarks. 



BASIC MODEL 

The basic model we consider is shown in Fig. |^(a) 
and comprises two identical ring BECs labeled j — 1,2 
which are in close proximity, and the whole system is 
rotating at an angular frequency ajR. The close prox- 
imity of the rings allows for spatially dependent tunnel- 
ing between them via mode overlap, meaning that the 
rings are coupled, allowing Josephson oscillations 
Each individual ring may be realized physically using a 
toroidal trap of high aspect ratio R = L/£o where L is 
the toroid circumference and £q the transverse oscillator 
length £q — yjfij muo , with loq the frequency of trans- 
verse oscillations, assumed to be harmonic. The trans- 
verse trap potential is assumed to be symmetric about 
an axis consisting of a circle on which the trap potential 
is minimum. The longitudinal (circumferential) motion 
on each ring can be described approximately by a ID 
coordinate Xj £ [— L/2,L/2] obtained by unfolding the 
ring and applying periodic boundary conditions, as illus- 
trated in Fig. 0(b). Then, at zero temperature the quan- 
tum dynamics of an atomic BEC moving on the paired 
rings may be described by the following coupled Gross- 
Pitaevskii equations in a reference frame rotating at lor: 
i i 
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where ipj^Xjt) is the macroscopic wave function for ring 
j = 1,2 with normalization condition 



dx{\M^,t)\'' + \Mx,t)\^) = N. 



(2) 



Here, N is the number of atoms of mass m, g — 
47r?i^a/(27r£Qm) = 2hLJoa > is the effective one- 
dimensional nonlinear coefficient describing repulsive 
many-body interactions, a being the s-wave scattering 
length, and il.{x) > 0, which is chosen real and positive, 
is the spatially dependent tunneling frequency between 
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FIG. 1: (a) The basic model we consider comprises two iden- 
tical ring BECs labeled j = 1,2 which are in close proximity 
and coupled via tunnelling, and (b) shows the unfolded rings 
to which periodic boundary conditions are applied. The rings 
come closest together at the origin a; = 0, where tunneling is 
represented by a dark oval. 



the two rings. In writing Eqs. (P we have taken ad- 
vantage of the fact that although the atoms in each ring 
are described by different coordinates a;j=i,2, they can 
nonetheless be described as moving on the same domain 
X e [—L/2, L/2] with the following caveats: First, the 
atoms on each ring do not cross-interact via mean-field 
effects, and are only coupled via the spatially dependent 
tunnehng. Second, inspection of Figs. 0(a) and (b) shows 
that atoms circulating from x = —L/2 — > L/2 along ring 
j = 1 are going counter-clockwise whereas atoms circu- 
lating from X = —L/2 — > L/2 along ring j = 2 are going 
clockwise. This means that although we write the equa- 
tions using a common spatial coordinate x G [—L/2, L/2], 
propagation in a given x direction corresponds to oppo- 
site senses of rotation for the different rings. This is why 
the rotation term proportional to luji in Eq. (|^) has a 
ring-dependent sign (—1)^ . 

With reference to Fig. |l|(a) we see, for example, that 
for an atom moving clockwise from a given reference 
point on ring j — I, then tunnelling over to ring j = 2 and 
moving counter-clockwise, and finally tunnelling back af- 
ter orbiting ring j — 2 to ring j = 1 to the original start- 
ing point, the atom crosses the tunnelling region twice. 
In this sense the coupled atomic rings are analogous to 
a symmetric SQUID j^], in which two superconducting 
rings are connected by a weak link, which has been em- 
ployed as a magnetometer [ p^ . The two ring system, 
then, is the simplest version of an Atomic-Phase Inter- 
ference Device, and we concentrate on the two-ring case 
in this paper to explore the basic properties of APHIDs. 

The tunneling frequency ^Imax — ^{x — 0) will be 
at its maximum at the point of closest approach of the 
rings, which we choose at a: = 0, and will decrease with 



separation, or equivalently as x varies away form the ori- 
gin. Typically, the tunneling frequency decays exponen- 
tially with ring separation. Thus, f2(a;) will typically be 
a bell-shaped function of x, and the spatial extent of 
the Josephson-coupling will be much less than the size 
of the ring L. Clearly, for smaller rings with tighter cur- 
vature, fl(x) will drop off faster away from the peak. In 
the limit Q{x) = 0, Eqs. (|^) reduce to the approximate 
one-dimensional form previously used to describe atomic 
BECs on a toroid. 

The conserved N-particle energy functional for the cou- 
pled Gross-Pitaevskii Eqs. (0) is 



E 



NhujQ 



dx[ 



2m 



5^1 



dx 



dx 



1^2 



hn{x){^ir2+rii^2)], (3) 



giving the energy per particle e = E/N. Since in this 
paper the transverse confinement energy huiQ is assumed 
the same for both rings and simply redefines the zero of 
energy, we hereafter drop this energy term for simplicity 
in notation. 



GROUND AND FIRST EXCITED STATES 

In this section we examine the properties of the ground 
and first excited states of a non-rotating (un = 0) pair 
of coupled ring BECs using a simple model to expose the 
main features. 



Zero-coupling limit 

It is useful in assessing the ground state properties to 
consider the non-coupled case with ^l{x) = 0. If all N 
atoms are homogeneously distributed on just one of the 
rings, with ipj — ^ N /L and — Oi then according 

to Eq. (^) the energy per particle is etrap — gn/2, where 
n = N/L is the linear atomic density. In contrast, when 
the atoms are equally split between the two rings, but still 
homogeneously distributed on each ring, [V'jl — \/ N /2L, 
and the energy per particle is 



gn 
4 ' 



(4) 



irrespective of the relative phase between the macro- 
scopic wave functions of the two rings. Energetically 
speaking then, in the absence of coupling the lowest- 
energy state is that in which the atoms are equally split 
between the rings as this minimizes the mean-field en- 
ergy. 
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Coupled solutions 

To proceed we now re-introduce the couphng and look 
for solutions where the atoms are equally split between 
the rings. In particular we consider solutions where the 
macroscopic wave functions of the two rings are in-phase 
(-I-) and out-of-phase (— ) by making the ansatz 



{±iy 

V2 



-^±*/V±(a:), J = 1,2 (5) 



with f±{x) the mode profiles on each ring and fj,± the 
corresponding chemical potentials. Then substituting in 
Eq. (pT) we obtain 



2m dx'^ 



W±?V±^TT'^{x)^±, (6) 



and J dx\(p±{x)\'^ ~ N . On general grounds, the out-of 
phase (— ) solution corresponds to the ground state. This 
can be seen from Eq. (|^) where the spatially dependent 
coupling Q.[x) > 0, which is typically bell-shaped, plays 
the role of a confining (de-confining) potential for the out- 
of-phase (in-phase) solution, thereby allowing for lower 
energy in comparison to the case without coupling. 

In the limit $1 = Eq. (™ also has the well-known 
dark soliton solution ||l3, 14, 1^, |l^, ^ |l^ on the infinite 



domain L ^ oo 



f±{x) = foix) — Vntanh ^ 



{X + Xq) 
V2: 



Xh 



(7) 



with /io = gn/2, where n is the linear density of the 
background (in the thermodynamic limit, where N —^ oo 
and L ^ oo, N/L — s- n remains non-zero). The healing 
length Xh is derived from the relation 



2mx^ 



gn 
2 ■ 



(8) 



The dark soliton solution represents a flat background 
density profile with a hole of width Xh << L located 
X — —xo, at which location a phase jump of tt also oc- 
curs as ipo goes through zero. In the thermodynamic limit 
the energy per particle associated with the dark soliton 
solution calculated using Eq. (||) is eo = ng/A — ei/2, 
that is, it is the same as that in Eq. (^) for a homoge- 
neous density on each ring without coupling. This arises 
because in the thermodynamic limit x/i/L — > 0, meaning 
that any energy increase due to the hole in the density 
makes a negligible effect on average; in other words, the 
hole in the density occupies a vanishingly small portion 
of the ring. 



Analytic approximation 

In general, numerical methods are required to solve Eq. 
|) for given parameters and tunneling profile il{x). In 



order to obtain insight into the ground-state properties, 
we employ a simple model 



fl{x) = il.max ■ d ■ S{x), 



(9) 



where flmax is the maximum tunneling frequency and 
d is the length of the tunneling region. This delta- 
function approximation will apply when d is much less 
than any other characteristic length scale of the problem, 
namely the ring length L and the healing length Xh- For 
the stationary coupled-ring solutions described by Eq. 
(^, where fl{x) plays the role of a single-particle poten- 
tial, the delta-function approximation yields a quantum- 
contact interaction |l^. Substituting Eq. (^) in (^ and 
integrating from a; = 0_ to a; = 0+ across the junction, 
we find that the action of the delta-function coupling is 
equivalent to a condition on the macroscopic wave func- 
tion derivative 



— ( g 

2m \ dx ^~°+ 



dx 



■d-'p±{0). 



(10) 

In the limit L » Xh » d we further impose the con- 
dition that ip±{x) is symmetric around a; = in order 
to satisfy the periodic ring boundary conditions, and we 
approximate 



ip±{x) 



m tanh 



{x + x±) 
V2: 



'Xh 



x>0. 



(11) 



With this approximation there is a cusp in ip±{x) at x=0, 
and the solution is extended to a; < by imposing reflec- 
tion symmetry around the origin. We can solve for the 
variables x± by substituting the approximate solution 
(|rT| ) in the boundary condition (|lO|), which yields 



■d = 



[1 - tanh^(a;±/V2a;,j)] 



\f2xh) tanh(a:±/A/2^ 



'Xh) 



(12) 



Since ^max > we find by inspection that the in-phase 
solutions correspond to x+ > and the out-of-phase so- 
lutions to x_ < 0. By introducing a dimensionless pa- 
rameter C = x+/ V2xh, with C > and C = x^/^/2xh, 
with C < 0, and using Eq. (||) for the healing length, we 
may write the above equation as 



■d 



9 



n [l-tanh2(C)] 
2n's tanh(lCl) 



(13) 



where Hs — mg/h? is a scaled density. Figure ^ shows a 
plot C versus the scaled tunneling frequency Mlmax ■ dj g 
for njns — 10^. Figure ^ shows examples of scaled den- 
sity profiles |v3±P/n for Ml^aax ■ d/ g = 'dbX — 0.5 (solid 
lines), hVl,nax ■ d/g = 6.3, C = 2 (dashed fines) and (a) 
the out-of-phase or ground-state solution, and (b) the 
in-phase solution. Density cusps in the solutions are 
evident, though we note that the ground-state density 
does not extend down to zero. The key features of the 
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Scaled tunneling frequency 

FIG. 2: Plot of versus h^lmax ■ d/g for n/us — 10*, with 
C = x+/V2xh,( > and ( = X-I^Xh,C, < 0. 

ground state are that as the scaled tunneling frequency 
ft^max ■ d/g is increased the depth of the density pro- 
file increases, the density at the origin going to zero as 
fi^max 'd/g —^ oo, and the width of the density hole also 
increases, approaching Xh as hflmax ■ d/g ^ oo. The in- 
phase solution is different in that it displays two density 
zeros and an on-axis maximum that is a cusp, shown in 
Fig. ^(b). Furthermore, inspection of the in-phase so- 
lution shows that its sign reverses through each density 
zero, and there are two sign reversals around each ring to 
ensure that the wave functions are single- valued. The in- 
phase solution therefore has a phase structure like a pair 
of dark solitons on each ring. For small Ml„iax'd/ g « 1, 
the density zeros are far apart (dashed line in Fig. ||(b) 
for Mlmax ■ d/g = 0.2, ^ = 2), but come together at the 
origin as fiVtmax 'd/g oo (solid line in Fig. ||(b) for 
fi^max ■ d/g = 3, C = 0.5). Thus, for both the in-phase 
and out-of-phase solutions the density vanishes at the 
origin as hflmax ' d/g ^ oo, and we have 

ip±{x) K y/ntSi-Dh [ /I ) . (14) 
\y/2xhJ 

A quantity of physical interest here is the energy per 
particle e± for the two solutions. Using the above ap- 
proximate solution in the energy functional (|^) we find 
in the thermodynamic limit 

e± = ^ ± ?jf^maa; • • ntanh^(ICI), (15) 

where the solution is again parameterized by C,. Note 
that in the limit of zero coupling ^max 0, the energies 
per particle of the two solutions become the same and 
equal to that of the equally split solution ei/2 = ng/4, as 




(a) 



(b) 



FIG. 3: Scaled density profiles \ip±\^ /n for n/us = 10*, 
hQrnax-d/g = 95, C = 0.5 (solid lines), TiQ,max-d/ g — 6.3, C = 2 
(dashed lines) and (a) the out-of-phase or ground state solu- 
tion, and (b) the in-phase solution. 

they should. Using Eq. ( p^ in ( p^ we obtain finally 

e± = ei/2 {\ ± y^tanh(|C|)[l - tanh2(C)]^ , (16) 

which is once again parameterized by C,. Figure |^ shows 
£±/ei/2 versus hQ.max ■ d/g for n/us = 10, the upper 
solid line corresponding to the in-phase (-I-) solution and 
the lower solid hue to the out-of-phase (— ) or ground 
state solution. For small values of the scaled tunnel- 
ing frequency hflmax 'd/g < 1 the energy per particle 
for the in-phase (out-of-phase) solution initially increases 
(decreases) away from ei/2 for zero-coupling, and this 
is expected physically. However, as the scaled tunnel- 
ing frequency is increased further the energy per particle 
for the in-phase (out-of-phase) solution reaches a turn- 
ing point at hflmax ■ d/g « 2, then decreases (increases). 
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FIG. 4: Scaled energy per particle e±/ei/2 versus TiQ.max -d/g 
for n/us — 10*, the upper solid line corresponding to the in- 
phase (+) solution and the lower solid line to the out-of-phase 
(— ) or ground-state solution. 



and both e± tend back to the zero-coupling value foi' 
fi^max 'd/g — *■ oo. The reason for this is that, as dis- 
cussed above, for both solutions the density tends to zero 
at the origin .t = where the junction is concentrated in 
the limit hflmax • d/g — > cx), so the Josephson coupling 
is rendered inoperative and the energy per particle tends 
to that for zero-coupling. 



V 



ng 



/3 = 



ng ng 

(20) 

These are the scaled equations used for our numerical 
study. We have solved the equations numerically using 
the split-step Fast-Fourier transform method [po| . 

To study the quantum dynamics of coupled-ring BECs 
we shall use an initial condition at r = where all N 
atoms are on one ring in a vortex state of winding num- 
ber p. This may be realized, for example, by condensing 
the atoms on one ring in the absence of the other, stir- 
ring the BEC to create the vortex and then turn- 
ing on the second ring. Sauer et al. have demon- 
strated a 2-cm diameter magnetic storage ring for laser- 
cooled, and Arnold and Riis (23j are working towards 
realizing a 10 cm diameter magnetically trapped toroidal 
BEC. One scheme for turning rings off and on is to use 
toroidal optical dipole traps formed by Laguerre- 
Gaussian beams piercing a two-dimensional BEC to cre- 
ate the rings j|, |2^, 26 , or alternatively using scanned 
laser beams to form the toroidal traps Cavity field 
enhancement may also be used to allow for large-radius 
toroidal traps ||2^. Regardless of experimental method, 
the initial condition we take is 



'^1(^,0) 



= 0. 



(21) 



Resonance conditions 



TIME-DEPENDENT SOLUTIONS 
Scaled equations 

For purposes of numerical simulations we introduce 
a simple Gaussian model for the spatially dependent 
Josephson-coupling 

ri(x) — ^max^ ' ^ ^max ' d • /(x), (1'^) 

with ^rnax thc maximum tunneling frequency and w << 
L the width of the coupling region around x = 0. 
We also introduce the normalized Gaussian f{x) — 
exp(— x^/w^)/ V TTw^ for which d — so that ef- 

fective parameters can be compared with the previous 
section. Then introducing the scaled variables 



(18) 



To proceed we examine the resonance conditions lead- 
ing to the initial exchange of atoms from ring 1 — > 2 using 
first-order perturbation theory. For the initial condition 
( pl| ) we choose the zeroth-order solution as that for v ~ Q 

(1)^^'^ r) = f,2^^pi^-m^'0p'-p^+i)r^ (22) 
Then writing the first-order solution for ring 2 in the form 



6« 



yields 



where 



q— — oc 



\aq{r)?-^,'^l _Bin(x,W2) 



pq 



Apq 



(23) 



(24) 



with n = N/L the mean density as before, we obtain 
with j = 1, 2, 



.9^ 

' dr 



2 9^2 



(19) 



Xpq = 2TT^/3{p^-q^)-iy{p + q) + l. 



(25) 



The vortex states q of ring 2 are therefore excited and 
generally exhibit small oscillations except at resonance 
where Xpq becomes small. The level of excitation of the 
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g*" vortex state is also dictated by the factor Tpq, but 
since we assume a narrow junction w/ L = A << 1, this 
factor allows for almost constant excitation Tr,, 



in the range q = p ± Sq with 



Sq 



1 

ttA 



>> 1. 



Consider first the case that the system is not rotating 
1/ = 0: Resonance occurs for that integer value of qr for 
which Xpq is equal to or closest to zero 



Qr 



P 



1 



27r2/3' 

the width of the resonance being 

1 
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(27) 



Aq: 



27r^p + qr)P' 



(28) 



When the width of the resonance is small Aq < 1 the 
initial vortex of index p in ring 2 will selectively cou- 
ple to vortices with mode indices qr satisfying Eq. ( |27| ) 
in ring 2, giving rise to relatively simple few mode dy- 
namics. In contrast, when Aq » 1 the initial vor- 
tex of index p in ring 2 will couple to a broad range 
of vortices with mode indices qr ± Aq in ring 2, giv- 
ing rise to multi-mode dynamics and complex behavior. 
In addition, for Josephson oscillations to occur the tun- 
neling energy per particle averaged over the ring length 
(1/L) J dx ■ hQ{x) = h^maxy^w/L should be for the 
same order as the mean- field energy per particle ng, or 



V 



ng 



V^A- 



(29) 



This gives an estimate of the scaled tunneling frequency 
r] to obtain Josephson oscillations. 



E 



o 



o 
50 



\ / / 

A 1 \ 

- \ \ 
\ 1 ' / 

\ ' 'i / 


\ '' / \ '' ' 

\ / \ \ \ 1 1 
\ ' * / ' ^ 


\ 1 \ 
\ ' \ 

V / 

- ;\ \ 
< \ / 


\ ; \ \ ; 11 
\ ' \ r 

' \ h ' \ 

/ \ / ■ / \ / " 


- / \ / \ 


\ / M / " 

' \ / ' ' \ / 
' \/ \/ 







1 

0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 
L 



(a) 



4 6 8 

Scaled time x 



10 




(b) 



4 6 
Scaled time x 



10 



FIG. 5: Fraction of atoms in each ring for A = 10 ,1] — 50, 
and (a) /3 = 1, Ag = 0.22, and (b) = 5.1 x 10"^ Aq = 1. 



Numerical results 

Here we present some examples of the dynamics of 
coupled ring BECs. For all the simulations we set 
p = 0, A — 10^2, and rj = 50. Consider first that the 
initial state corresponds to the ground state {p = 0) of 
ring 1. From Eqs. ( p7| ) and (|2^) we obtain 

^r-^^=^,, (30) 

that is, the width of the resonance Aq is equal to the 
resonant value q = qr- Figure ^(a) shows the fraction of 
atoms in each ring for (3 — 1 for which qr — Aq = 0.22, 
and complete Josephson oscillations between the two 
rings are evident. In this case the density profiles in 
the two rings are largely flat as resonant coupling occurs 
between p = 0, ~ 0. In contrast, for /3 = 5.1 x 10^^ as 
shown in Fig. |^(b) for which qr = Aq — 1, the Joseph- 
son oscillations are now incomplete. Physically, there are 



multiple modes involved in ring j = 2 with g = 0, ±1, ±2, 
and the resulting multi-mode dynamics is what frustrates 
the Josephson oscillations for Aq > 1. The multi-mode 
dynamics manifests itself as spatial density modulations 
in the two rings as shown in Fig. |^(a) for the same pa- 
rameters as in Fig. |^(b) and r = 10. For even lower 
density /3 = 5.1 x 10""* for which qr — Aq = 10, the 
Josephson oscillations are all but extinguished, and the 
spatial density proflles in rings j — 1,2 are shown in Figs, 
^(b) for r = 10. Clearly, the multi-mode nature of the 
solution allows the coupling due to tunneling to concen- 
trate around the coupling region, hence reducing the net 
fraction of atoms transferred between the rings. 

Some estimates of parameters are in order. Using 
g — 2huJoa gives ng — 2%uiQN{a/ L), and Vlmax — 
2TjLUoN{a/L). Then for = 27r x 10^ rad-s-\7V = 
10'^, L ~ 1 cm, a = 5 nm, we find flmax — 2Tr x 5 rads^^, 
and T is time in units of fi/ng — 1.6 s, so the Josephson 
oscillations in Fig. ||(a) occur on a time scale of seconds. 
Setting m = 10"^^ kg we obtain = mg/Ti' k. 63 cm~^. 
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FIG. 6: (a) Spatial density in ring j = 1 (solid line) and 
j = 2 (dashed line) for the same parameters as Fig. ^(b) with 
T = 10, f3 = 5.1 X 10~^, and Aq = 1; (b) spatial density 
in ring j — 1 (solid line) and j = 2 (dashed line) for /3 — 



5.1 X 10"*, Ag = 10 for r = 10. 



and for n = N/L = 10^ cm^\ p = 1.6 x 10"^ It is im- 
portant that n/us > 1 to ensure that the one-dimensional 
gas acts as a BEC as opposed to a Tonks gas |3^ . The 
parameter f3 = (n/ns)/N'^ is proportional to I/ujq and 
1 /N so we can increase /3 by decreasing either the num- 
ber of atoms and/or the transverse oscillator frequency 
with respect to the above values. 



Effects of rotation 

An interesting feature of the two ring APHID is that 
the condition Xpq = 27r^/3(p^ — q^) — v{p + g) -I- 1 — > for 
resonant coupling between the rings is dependent on the 
scaled rotation rate u — hLUn/ng. In particular we find 
for p = 
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FIG. 7; Percentage of atoms in ring j = 2 as a function of 
scaled rotation rate u — hcjR/ng. 



(31) 



This implies that for scaled rotation rates > y^Sn^fi 
the rotation of the entire APHID will affect the coupling. 
Consider then a case where without rotation Aq >> 1 so 
that the Josephson-oscillations are all but extinguished 
and the atoms remain on ring 1. Then as the scaled 
rotation rate f is increased from zero, inspection shows 
that one solution qr in Eq. (^) moves towards reso- 
nance while the other moves further away. Therefore, 
starting from a detuned case with minimal coupling, in- 
creased rotation leads to increased coupling which can 
then be detected via the number of atoms on ring 2 at 
a fixed detection time. Figure |^ shows the percentage of 
the atoms in ring 2 versus scaled rotation rate v at time 
T = 10 and A = 0.01, k = 50, /3 = 5.1 x lO"", for which 
q,. — Aq — 10, and the effect of rotation dependent cou- 
pling between the rings is clearly exhibited. Some points 
are worth making here: First, the rotation causes the 
number of atoms in ring 2 to change by about 10% of 
the total number of atoms, so experimentally it will be 
necessary to control the initial number of atoms on ring 
1 to better than this percentage. Furthermore, it would 
be a challenge to detect the small number of atoms in 
ring 2. Second, for our particular example with p — 
the number of atoms in ring 2 is sensitive to the mag- 
nitude but not the sign of the rotation, but this can be 
changed by having p ^ in which case the the coupling 
becomes sensitive to the sign of h'. Third, the sensitivity 
of the atom number to rotation rate increases with the 
observation time t chosen, remembering that we are in a 
far-off-resonant situation so coupling happens slowly. Fi- 
nally, the number of atoms in ring 2 is not necessarily a 
monotonic function of the rotation rate, as seen from Fig. 
0, which will limit the range of rotation rates that can 
be uniquely measured. Nonetheless, we feel this is an in- 
teresting phenomena which may have utility for rotation 
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sensing with further development. 

To gain some sense of the sensitivity of this scheme we 
use the same parameters as the previous section for which 
r is time in units of h/ng = 1.6 s. Then a value of = 
0.01 corresponds to a rotation rate ujr = 2i'U!oN{a/ L) — 
2tt X 10^"^ rads^^ which is one hundred times higher than 
the Earth's rotation rate at the poles. However, if we are 
willing to reduce the transverse oscillator frequency to 
ujQ — 2tt X 1 rads^^, then ly — 0.01 corresponds to the 
Earth's rotation rate, but then time is in units of 160 
s in the figures! We are currently working on schemes 
involving multiple-ring APHIDs to enhance the rotation 
sensitivity. 



SUMMARY AND CONCLUSIONS 

In summary, we have presented a theoretical investiga- 
tion of a pair of ring BECs coupled by tunneling as the 
simplest example of a potential Atomic Phase Interfer- 
ence Device. We have shown that the two-ring APHID 
has interesting ground-state properties, with density pro- 
files reminiscent of dark soliton states around the point of 
contact of the rings. Furthermore, we have demonstrated 
that Josephson oscillations between the two rings can oc- 
cur, and that these oscillations are sensitive to the state of 
rotation of the APHID. In particular, if all the atoms are 
prepared on one ring, then the number of atoms trans- 
ferred to the second ring in a given time span is a measure 
of the rotation rate of the APHID. Although the two-ring 
APHID was found to be not very rotation sensitive, we 
believe APHIDs are worthy of further study as multi- 
ring APHIDs will display enhanced sensitivity to the rel- 
ative phase between the rings hence potentially leading 
to increased rotation sensitivity. We shall be exploring 
multi-ring APHIDs in future research. 
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